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Abstract 

In this paper, we present an approach for dynamic exploration and mapping of unknown 
environments using a swarm of biobotic sensing agents, with a stochastic natural motion model 
and a leading agent (e.g., an unmanned aerial vehicle). The proposed robust mapping technique 
constructs a topological map of the environment using only encounter information from the 
swarm. A sliding window strategy is adopted in conjunction with a topological mapping strategy 
based on local interactions among the swarm in a coordinate-free fashion to obtain local maps 
of the environment. These maps are then merged into a global topological map which can 
be visualized using a graphical representation that integrates geometric as well as topological 
feature of the environment. Localized robust topological features are extracted using tools from 
topological data analysis. Simulation results have been presented to illustrate and verify the 
correctness of our dynamic mapping algorithm. 


1 Introduction 

Sensor networks with their broad application in mapping and navigation , habitat monitoring , 
exploration, and search and rescue [^, have attracted a lot of attention in recent decades. Mobile 
sensor networks offer the flexibility to adapt with dynamic environments. In some applications, 
however, the amount of information that could be sensed or transferred by the agents is limited. 
This motivates the design of distributed systems composed of simple agents with minimal sensing 
requirements [^. 

As a motivating example, imagine a scenario in which an earthquake has left several individuals 
trapped under the ruins of collapsed buildings. Sending teams of human rescuers to search and 
extract survivors may put the rescuers and the survivors in danger. Instead, we can choose to 
send multiple autonomous robotic agents to accomplish this search and rescue operation [^. Since 
large robots may not be able to explore all desired locations due to size and safety limitations, 
a team of small autonomous agents can be considered (e.g., biobotic insects II1I!I or biologically 
inspired milli-robots [^; see Fig. [^for illustration). The first requirement for such system would 
be to explore and build a map of the environment and then localize anyone in need of assistance. 
This is but one scenario in which a swarm of small robotic agents can be used in unstructured 
emergency-response situations for mapping and source-localization purposes. 
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Figure 1: A physical environment with a swarm of biobotic (cyborg insect) agents on the ground with 
their corresponding local sensing neighborhoods and an aerial leading agent. 


The mapping task in such situations becomes extremely challenging due to hardware limitations 
and the unstructured nature of the environment. Power and computational resource constraints 
prohibit us from using continuous control schemes for the agent’s locomotion and from using on¬ 
board imaging techniques for their localization. Furthermore, since some locations can be indoor 
or even underground in cluttered environments, signal propagation based localization (e.g., GPS, 
or computing signal strength or time of flight) may be unreliable, and odometry information might 
include a high amount of uncertainty due to irregular conditions of the terrain. Therefore, tradi¬ 
tional mapping and exploration techniques such as SLAM would not perform well under these 
adverse conditions. 

Since obtaining an accurate metric map of the environment may not be possible in these scenar¬ 
ios, we are interested to construct topological maps of unknown environments using biobotic mobile 
sensor networks under the constraint of limited sensing information. Our approach to mapping in 
such scenarios is to extract a sketch of the environment instead of a fully-detailed map. This sketch 
is composed of robust topological information that is obtained from a minimal amount of sensing. 
Instead of providing continuous control feedback to the agents, we explore how stochastic motion 
models and weak encounter information can be exploited for mapping. These strategies for motion 
and sensing accommodate for the hardware limitations of the platforms under consideration. We 
employ tools from algebraic topology to extract spatial information of the environment based on 
neighbor to neighbor interactions among the agents with no need for localization data. This in¬ 
formation is used to build a map of persistent topological features of the environment, which will 
contain only robust features. 

In this paper, we extend the approach introduced in 


10 , 11 


where a notion of encounter 
metric is built upon encounters among the agents, and uses topological data analysis (TDA) [12| 


to extract features such as connected components and holes in the environment. This metric was 
improved in 11 to obtain a more accurate point cloud representation of the unknown environment. 


by using a density based subsampling procedure as well as a classification technique to provide a 
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Figure 2: Illustration of the physical environment to be explored with the swarm of agents performing the 
exploration and mapping, where A is the area of interest and 'D{t) denotes the covered domain by the aerial 
agent flying with the velocity of u{t) at time t. 
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robust quantitative representation of features in the point cloud data. Although the aforementioned 
approach works quite well in simple scenarios, its performance in accurately classifying features 
degrades for more complex environments, with growing number or size of persistent features in the 
environment [^. Additionally, the complexity of the algorithm grows quadratically with respect 
to the number of observations required to properly estimate a map of the scene. This makes the 
procedure computationally too expensive for larger and more complex environments. 

Therefore, in this paper, to avoid increasing the computation cost and to obtain more accurate 
and robust maps for larger environments with more complex features, we propose a dynamic ex¬ 
ploration and mapping strategy by addition of a leading agent (see Fig. and [^. This approach 
can be described as a two level procedure: Globally, at a higher level, the aerial agent performs 
a dynamic sweeping coverage of the area of interest by leading the center of the swarm to move 
towards a desired direction. Locally, the swarm of boibotic agents explore their assigned coverage 
domain. Information gathered over smaller time windows are exploited to build local maps of the 
domains covered and then merged into a global map (Fig. [^. 

The remainder of the paper is organized as follows: Section overviews the related work to 
this study. A concise background on the tools from topological data analysis used in this work is 
provided in section section describes the problem under study, including mobility and sensing 
characteristics of the ground biobotic agents as well as the leading agent; in sectionwe present our 
methodology for dynamic exploration and topological mapping with minimal sensing information; 
section 1^ validates our approach using numerical simulations; hnally, conclusions and future work 
are discussed in section [7l 


2 Related Work 


The standard localization and mapping techniques, originated in the SLAM algorithm [^, estimate 
the position of moving agents based on a set of observations from landmarks in the environment as 
well as sensing information from the moving agents (e.g. odometry, IMU, or vision) to estimate the 
agents’ motion 13 . As an example, a lot of recent studies on localization and mapping, has focused 
on visual or RGB-D SLAM |14l|15|. The primary assumption in these approaches is, however, that 


the exploring agents are equipped with advanced sensing devices such as depth sensors or camera. 
In our problem, on the other hand, we consider a case where such sensing information is unavailable 
(e.g. biobotic or micro-robots), which will result in the failure of standard methods. 

Methods from computational topology, on the other hand, can provide tools to extract topo¬ 
logical features from data sets without requiring coordinate information. This makes them more 
suitable for scenarios in which weak or no localization is provided. Topological frameworks, which 
apply tools from algebraic topology, have been recently used for a variety of applications including 
sensor networks and robotics. Particularly, these methods have been used for applications such as 
coverage and hole detection in sensor networks 


16 -20 , motion planning 21 , localization p2 -24 


and camera network coverage [25] . 

Walker [26j employed persistent homology to compute topological invariants from encounter 
data of the mobile nodes in Mobile Ad-Hoc networks in order to infer global information regarding 
the topology of a physical environment, but the nodes are assumed to follow a simple mobility 
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model on a graph. A topological version of localization via signals of opportunity has been studied 
in 24 . However, this approach is mostly focused on localization problem although it could also 
detect coarse features of the environment in the context of mapping. Furthermore, it requires the 
existence of transmitters in the environment that provide such signals of opportunity. Topological 
localization and mapping has been also considered lately in |22| and |23| . Although these are akin 
in the context of using topology for representation of physical structures in the environment with 
mobile agents, they differ with our work in the sense that they assume the existence of an IMU or 
a motion model for prediction of the agents’ states, as well as fixed landmarks on the boundaries 
of structures (e.g. buildings) to be used for observation. In this paper, however, the mapping is 
intended to be performed using local encounters among the agents and not relying on external 
resources (signals of opportunity or fixed landmarks), as they might get out of order as a result of 
disaster scenario conditions, nor on odometry information. 

In 1^,27 , topological characterization of coverage and hole detection in sensor networks using 
only proximity information of the nodes within a neighborhood was examined. However, these 
studies focused on stationary networks, and mainly concerned about the coverage holes in the 
sensing domain of the network rather than the characterization of the physical environment itself. 

There exists also a notion of topological SLAM in the literature (e.g. 


28,29 ), which differs 


from our work in the sense that the context of topology has been considered as a combinatorial 
representation of the free space of the environment. 


3 Background 

3.1 Toppological Data Analysis 

A brief introduction of some of the basic concepts in computational topology and persistent ho¬ 
mology used for data analysis in this paper is presented here. A comprehensive review of the topic 
can be found in 30 . 

Topological data analysis (TDA), introduced in [^, is a new field of study which employs 
tools from persistent homology theory (Ml to obtain a qualitative description of the topological 


attributes and visualization of data sets sampled from high dimensional point clouds. 

A point cloud can be thought of as finite samples taken from a density map which may in¬ 
clude noise. TDA represents the prominence of features in the point cloud in terms of a compact 
representation of the multi-scale topological structure called persistence diagrams (^. It reduces 
the dimensions of data by construction of a filtration of combinatorial objects, which can represent 
geometrical and topological features of the data set at specific scales, and the nodes at this complex 
can be considered as corresponding clusters in the data. 

One of the well-known techniques widely used in topological data analysis is persistent homology, 
which deals with the way that objects are connected. Topological structures of a space M are 
summarized as a compact representation in the form of so-called Betti numbers, which are ranks 
of topological invariants, called homology groups. The n-th Betti number, /?„ measures the number 
of n-dimensional cycles in the space (e.g. /3o is the number of connected components and /3i is the 
number of holes in the complex). The space M usually is not directly accessible but a sampled 
version of it, X, can be used for computations. This sample is represented as a point cloud, a finite 
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Figure 3: Topological persistence: (a) An example of a topological space M with a sampled point cloud 
X, (b) a filtration of simplicial complexes over X, (c) the corresponding persistence diagram with dgmo(A) 
in blue and dgmi(A) in red, highlighting persistence of connected components and holes over the scale, 
respectively. 


set of points equipped with a metric, which can be defined by pairwise distances between the points 
(distance function). The distance function can be defined based on the coordinates of the points in 
their embedding manifold, or can be constructed by specifying pairwise distances between points 
independently from their (possibly unknown) coordinates. 

A standard method to analyze the topological structure of a point cloud is to map it into 
combinatorial objects called simplicial complexes. One way to build these complexes is to select 
a scale e, place balls of radius e on each vertex, and construct simplices based on their pairwise 
distance relative to e. A computationally efficient complex, called the Vietoris-Rips complex 32 , 
consists of simplices for which the distances between each pair of its vertices are at most e. A 
sequence of complexes, called a filtration Ai(e), then can be obtained by increasing e over a range 
of interest, with the property that if f < s then X{t) C -T(s). Persistent homology, computes 
the values e for which the classes of topological features appear (6(j) and disappear (dj^) during 
filtration, referred to as the birth and death values of the i-th class of features in dimension n. This 
information is encoded into persistence intervals [6^, dl2 or as a multiset of points {bl^, d\), called a 
persistence diagram, dgm„(X). Algorithms for computation of persistent homology can be found 


in 12,33 


An example of a topological space M together with sampled data X and corresponding filtration 
over e is shown in Fig. The persistent diagram infers the existence of one connected component 
and one persistent hole. 


4 Problem Statement 

Consider a collection of mobile sensing agents moving stochastically on a bounded domain of interest 
C M^. Each ground agent is distinguished by its unique ID, and its motion dynamics in the 
bounded space mimic the movements of a biobotic insect. In addition to the swarm of ground 
agents, there exists a leading agent (e.g. a quad-copter or a biobotic moth), which plays the role of 
a swarm leader in order to herd the ground agents to explore all parts of the desired environment 
over time. It serves as an access point to establish a local network for communication with the 
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swarm in order to send control commands and receive required information (Fig. [^. 

The ground agents are provided with weak localization information, i.e. neighboring agent 
identification within a specific range are the only information used, and no coordinate information 
or other information is provided to them. For example, odometry and inertial measurements 
are considered to be too unreliable due to the uneven and unstructured terrains present in our 
application scenario. However, they can communicate wirelessly to the leading agent within its 
limited communication range. Alternatively, due to weak transmitters on the ground agents, we 
could just consider a one-way communication from the aerial agent to the swarm. 

Given this minimal sensing scenario, we aim to build a sketch of the unknown scene, which 
includes information about topological and geometrical features of the environment A. We refer 
to such sketch as the topological map. Fig. [^illustrates the idea of such topological mapping for a 
square shaped physical environment including two obstacles. Our description incorporates a point 
cloud representation obtained from local interactions between agents (Fig. j^b)), a list of robust 
features captured by topological persistence diagrams (Fig. Qc)), and a physical sketch of the 
environment (Fig. j^d)). 

The previous example demonstrates the approach proposed by the authors in [^, in which 
stochastic agents completely explore a bounded domain in order to construct a topological map. 
These results will be the building blocks for the approach proposed in this paper, in which an 
external agent is used to lead the swarm in order to allow for mapping and exploration of large 
and more complex environments. We assume the existence of a leading agent (e.g., an aerial 
vehicle) which sweeps the entire domain. The stochastic agents are assumed to stay within a 
region T){t) C A centered around the leading agent. As the swarm of agents travel through the 
environment, a dynamic map of the scene is built by having the agents share information about 
their interactions. In order to provide enough encounters between agents to be gathered, the leading 
agent must travel at an appropriate speed. Furthermore, we require that the dynamic coverage of 
the stochastic swarm to cover the domain of interest. That is. 


where [0, T] is the observation time interval. 


( 1 ) 


4.1 Motion Model of the Biobotic Swarm 

The mobility model of the ground mobile sensors is adopted from the probabilistic movement 
model of cockroaches, described in [34] ,which can be mainly described by their individual and 
group behaviors. 

The individual behavior of cockroaches can be described mainly by two behaviors, namely 
diffusive random walk (RW) and wall following (WF). In this paper, for the sake of simplicity, 
we skip the group behavior and wall following and only consider agents who stochastically move 
according to the RW model. This is consistent with the behavior of cockroaches 
exploited as a natural mode of operation for biobotic systems . 

The random walk is modeled as piecewise linear movements with fixed orientation, interrupted 
by isotropic changes in direction, and constant average velocity of Vm- The length of the line seg- 


34 , which can be 
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Figure 4: Topological Mapping: (a) physical environment, (b) estimated point cloud from encounter infor¬ 
mation, (c) persistence diagram highlighting connected components in blue and holes in red dots (significant 
features can be distinguished from noise using an appropriate threshold depicted as a red line), and (d) 
sketch of the environment. 


ments are selected randomly from an exponential distribution with characteristic length A. Changes 
in direction are also triggered by collision detection with other agents or obstacles in the scene. 

Moreover, it is known that during their RW motion, the agents probabilistically stop for some 
period of time and then continue their movement [^. The biobotic insects are equipped with a 
controlled motion mode aside with their natural motion mode. In the controlled motion mode, the 
insects are capable of being controlled wirelessly to perform left and right turns or to be commanded 
to start or stop their motion. Consistent with this controlled mode, we assnme that, aside from 
random walk, the agents can be in a static mode of operation in which they do not move but sensing 
and communication are still possible. The only control that we assume on their motion is switching 
between static and random walk modes. 

Note that the average velocity of the swarm is determined by the average velocity of the flying 
leader, u. In order for the biobots to fully explore and sample the current coverage domain of the 
leader with high enough density, the average velocity of the random walk motion of the swarm 
agents needs to be high enough compared to the swarm velocity. 

4.2 Sensing and Communication 

The sensing and communication model for the agents is inspired by limited sensing capabilities of 
the biobotic insects, which incorporate a wireless transmitter and receiver provided by a system- 
on-chip based ZigBee enabled wireless backpack system inserted into their bodies . Each agent is 
distinguished by its unique ID, and assnmed to have a limited sensing capability, with a detection 
radius r^, which defines the region in which it can identify and communicate with other agents. It 
is assumed that the swarm can only detect signal strength from the leading agent, which can be 
used to maintain the agents within a bounded domain D{t) centered around the leading agent at 
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Figure 5: Dynamic topological mapping (a) A block diagram illustrating several components of mapping, 
grouped into three major tasks: exploration, local mapping, and global merging, and (b) components of 
local map construction over each time window. 


time t. 

Moreover, each agent can detect the presence of an obstacle and avoid collision by changing 
direction within its collision detection region. The agents are able to record their encounters with 
each other in accordance with the corresponding times as encounter events. These encounters are 
in fact the main piece of information in our approach in the construction of a metric on the samples 
of the physical space for the purpose of mapping. An encounter event takes place between two 
agents if their corresponding center positions fall within a distance of r^. An encounter Tj can be 
recorded by an agent as a tuple: 

£i = [ioi, tii, ID^], (2) 

where [toiTii] denotes the time interval of the encounter, and and ID? represent the IDs of 
the two encountering agents. We define T{£i) = [toijiii] as the i-th event time interval, and 
ID{£i) = {ID^IDf} as its ID set . At each time step, we consider a subset of agents in the 
corresponding coverage domain P(t) to be in the static mode, and denote by S{t) the set of static 
nodes at time t. Furthermore, the nodes are able to report their status as being in a random walk 
or static state. 


5 Dynamic Topological Mapping 


As described in the introduction, our objective is to construct a map of the environment in scenarios 
with no reliable odometry, signal strength or visual features (due to an unstructured and harsh 
terrain), and with limiting power constraints. For this reason, we rely only on a minimal sensing 
strategy that makes use of encounter information between stochastic agents (i.e., keeping track of 
interactions between agents that are in close proximity of each other). 

In this section, we represent our methodology to extract topological and geometrical features 
from unknown environments. An overview of the approach is shown in Fig. The basic method¬ 
ology makes use of the topological mapping approached introduced by the authors in 


10 


as a 


building block. As the leading agent sweeps the environment, observations over a set of windows of 
time are used to construct local topological maps. Information from observations in the intersection 
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of these windows are used to determine connections between these maps in order to stitch them 
together and build a global topological map. 

We define At as the total amount of overlap between two time windows. This overlap needs to 
be sufficient in order to guarantee that enough static points and encounter events are found there 
to be able to properly stitch the topological maps (see section 5.2 for more details). The time 
windows from which observations will be used for the local mapping are defined as 


Wi = [ti-i - At/2, ti + At/2] = (3) 

where {tj} are evenly spaced, N is the total number of windows, and IJi = [0)^] (see figure]^ 
for an illustration). 

As mentioned before, D{t) is the domain of interest around the leading agent located at x{t) in 
which the ground agents will be exploring the domain (e.g., D{t) can be a disk of radius r centered 
at x{t)). Hence, we can define the i-th domain for local mapping as 

A = U D{t), (4) 

teWi 

which is the region in the space to be mapped locally. In order to guaranteed full coverage of the 
environment A, we must have A C 


5.1 Local Mapping 

This task is accomplished by following the methodology introduced in [^. The ground agents 
perform the exploration and mapping by simultaneous random dispersion and leader following. In 
order to create fixed landmarks in the environment, a subset of agents are selected as static. The 
process is executed by: (1) recording encounter events between agents, (2) constructing a metric 
based on these events, and (3) extracting robust topological features. 


5.1.1 Exploration and Event Collection 


The exploration initiates with a dispersion stage, where all nodes start their motion from an 
arbitrary initial configuration in a moving status until dispersed throughout the environment, while 
maintaining their connectivity with the aerial access point. This is followed by a landmark selection 
stage, in which a small percentage of the nodes are selected as static landmarks by modifying 
a standard MaxMin landmark selection approach 35 . Given the adjacency graph due to the 


communication range, we use its hop distance to select landmarks iteratively by maximizing its 
distance from the previous set of landmarks. The first agent is selected randomly. If the graph is not 
connected we assume that landmark agents are selected randomly from each cluster with probability 
equal to the percentage of agents in each cluster. We have shown in |36| that incorporating a small 
percentage of static nodes in the network, as landmarks will improve our point cloud estimation of 
the environment to a great extent. Once the static landmarks are selected, the rest of the agents 
continue their random walk motion and record their encounters with other agents. 

In order to incorporate the necessary aspects of leader following, we require the ground agents 
to remain within T)[t) for all times. We assume that the agents are capable of detecting when 
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they leave this domain of interest, and can change their random walk strategy in order to estimate 
a local signal strength gradient and move in that direction. Static agents that leave 'D[t) would 
immediately switch to random walk behavior, and new landmark agents would be selected. Our 
method is robust in the sense that it tolerates having some agents getting separated from the 
swarm. 


5.1.2 Geometric Reconstruction 

As our approach is a coordinate free procedure, we need to hrst build an estimated metric to 
represent the map of the local region as a point cloud. Our metric, which we refer to as the 
encounter metric, is constructed based on a set of encounter events among the agents. 

To build a distance metric on the set of encounter events, we construct an undirected weighted 
graph Q with vertices corresponding to the events £i, denoted as the encounter graph. For any two 
vertices £i and £j, we include an edge Cij with weight 

Wij = inf {\ti - t 2 \ : ti e T{£i), t 2 e T{£j))} (5) 

if ID[£i) n ID{£j) 7 ^ 0. We can furthermore improve our graphical model by incorporating the 
fact that landmark nodes are static. Let S be the set of indices of all agents who has been static 
for some time period. For each s G 5 let 

r(5) = {ri(s),r2(s),...} (6) 

denote the set of all time intervals r*(s) = [tl’^, tl’'^] for which agent s has been in a static mode. 
Then, two events £i and £j with 

ID{£i) n ID{£j) = s£S (7) 

had to occur geometrically at nearby locations if 

3k,s.t. T{£i)nT{£j)nT'^{s) (8) 


Therefore, the corresponding weight Wij is set to 0 in this case. 

Note that the elements of the graph Q, represent estimations of pairwise distances between 
encounter events, i.e. G{i,j) = d£{£i,£j). To proceed with the geometric point cloud reconstruction, 
we follow the Isometric mapping (Isomap) technique [37] in manifold learning 38 . The main idea 


is to estimate Geodesic distances on a manifold by the shortest path distance on the neighborhood 
graph of the samples taken from the manifold. 

Therefore, we build our encounter metric on G, denoted by Vg, as [T’gJij = dg{£i,£j) where 
dg{£i,£j) represents the length of the shortest path between nodes £i,£j in G, which can be ob¬ 


tained from Dijkstra’s or Floyd’s shortest path algorithms |39| . Finally, we apply classical Multi¬ 
dimensional Scaling (MDS) to the matrix Dg in order to find an embedding of the point cloud 
data in 3-dimensional Euclidean space. 
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5.1.3 Topological Feature Extraction 


The set of encounter nodes together with the corresponding estimated metric defines a point cloud, 
the encounter point cloud, which is processed to construct a filtration of simplicial complexes, 
denoted as encounter complexes. This filtration is obtained by computing the Vietoris-Rips com¬ 
plexes [32j from the encounters as a function of the distance parameter e (as defined in section]^. 
The construction of the encounter complex is done in a computationally efficient and robust way 
by incorporating strategies as described in [IT, 36 


We use the tools introduced in section in order to extract persistence diagrams from a filtration 
of the subsampled set of points. The filtration is based on computing the Rips complexes [^ and 
we use the Dionysus C-|—|- library [^ for computation of persistent homology. The diagrams 
capture the persistence and robustness of topological features in the environment. We restrict our 
computations to the dimensions 0 and 1 of persistence as the higher dimension are not applicable 
in our experiments. Robust classification of topological features in the environment is performed 
by using statistical learning techniques as described in [TI]. 


5.2 Global Mapping 

So far we have constructed a set individual local maps {M.i} for the covered domains {Vi}. The 
objective of this section is to provide a methodology for gluing these maps together to obtain a 
global one, M. for the area of interest A. Figure illustrates our proposed approach. The idea is 
to use common static nodes between different map components to define a connectivity measure 
and connect them based on the defined measure. To preserve topological features that cannot be 
described on a single local map (e.g. a very long building that will not fit within one window), we 
construct point clouds on the intersection of the consecutive map components (M* and Mj+i) and 
perform a clustering analysis to robustly identify the number of connections between the different 
map components. 

We define the subset of static nodes joining Adi to Mi+i as: 

= {s G 5 : 3k s.t. <ti — At and > U + At}. (9) 

This is the set of static nodes whose appearance time precede the starting time of ITj+i and their 
disappearance time is later than the end time of Wi. Furthermore, we define the set of encounter 
associated with these static nodes as 

= {£i ■ ID{£i) n A 0}- (10) 

An encounter metric is constructed on the set Ssa+i to approximate a metric on the set Si^i+i. 
We consider a point cloud, consisting of the set T 5 . together with the k-nearest neighbors of 
each of its elements as a connecting point cloud between the maps Aii and Afj+i. We refer to this 
point cloud as the inter-domain point cloud and the equivalent metric as the correspondence metric 
between two local maps Aii and Mj+i. 

Note that Si^i+i correspond to connections between the two maps, but some of these connections 
may correspond to the same connected component within the domains Di and Dj+i. In order to 
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identify which connections are equivalent, we perform a clustering analysis by using the encounter 
metric constructed on the set Ssi j+i • 

As we do not know the number of clusters a priory, we employ the single linkage clustering as one 
of the most famous agglomerative hierarchical clustering techniques [^. Hierarchical techniques 
provide a similarity dendrogram which can be used to classify data into clusters by applying a 
cut-off threshold c on the similarity values. This cut-off is selected experimentally. 

6 Simulation Results 

In this section, we present the simulation results for the verification of our dynamic topological 
mapping approach. We consider basic geometries in our analysis in order to provide the building 
blocks for a structured mapping of more complex environments. In our simulations, we consider 
three different scenarios which differ in the physical layout of the area of the interest, as shown 
in the first rows of the Figures [^ - [^ The idea behind these scenarios is to verify the successful 
implementation of the proposed approach when there are large features (e.g. long structured 
obstacles) that cannot fit within a local coverage domain, like the one in Fig. or scenarios like 
in Fig. in which although the individual features are not too big, but need to be recognized and 
merged properly. 

The aerial leader is considered to move with a constant linear velocity of Vx = 0.005 m/s in 
the X— direction, while the average velocity of the ground agents in their random walk model 
is considered to be Vm = O.lm/s. The coverage window of the aerial vehicle, for simplicity, is 
considered to be a rectangular arena with the length of I = 3m, and the same width as the 
environment, where the dimensions of the environments are set as 10m x 30m. In each figure, the 
hrst row represents the physical environment together with the initial condition of the agents (gray 
circles) at t = 0 as well as the location of all static nodes (landmarks) over all time periods (shown 
as pink dots). For each scenario, a total number of 100 agents are considered in the network, among 
which 10% are selected as landmarks after some initial dispersion based on the maxmin landmark 
selection algorithm. The detection radius of each agent for encounter and boundary recognition is 
adopted to be = 1cm. 

We divide the whole run time of scenarios 1, 2, and 3 into 4, 6, and 5 windows, respectively, 
and implement our proposed local map construction on gathered encounter information over each 
window. Note that such selection have been made for sake of better visualization of the features, 
but in general we could select the same number of windows for all scenarios. The point clouds 
constructed from the estimation of local encounter metrics for different windows are presented 
in the third row of each figure (in pink color), in accordance with their persistence diagrams 
in the second rows. On each persistence diagram, 0-dimensional and 1-dimensional topological 
features (connected components and cycles) are plotted as blue and red dots, respectively, and 
the corresponding thresholds for the persistence classihers are sketched as dashed lines. It can 
be observed the point clouds representing local maps represent the correct topological features 
of the corresponding local domains covered over each window. This can also be verified by the 
topological features classified as robust ones in the corresponding persistence diagrams. Moreover, 
the estimated point clouds also contain some geometric information, which could be exploited to 
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improve the mapping accuracy in future work. 

In the next row of each the figure, the inter-domain point clouds obtained from the correspon¬ 
dence metric between local maps are plotted followed by the dendrograms representing the results 
of the hierarchical clustering on each of them. To specify the number of clusters for each point 
cloud, a clustering threshold value of 50 has been selected for all scenarios. In all cases, except 
the last three point clouds in scenario 2, they are classihed as one cluster. On the other hand, 
in scenario 2, as a result of the large obstacle, local maps A ^4 and and their corresponding 
inter-domain point clouds are all divided into two clusters. Finally, from the collection of local 
maps, and obtained correspondence metrics, we stitch the consecutive pieces of the maps together 
in order to find a global topological representation of the environments. The final result of such 
maps for each environment are sketched in the last rows of the Figures [^-|^ 


7 Conclusion and Future Work 


In this work, we presented a dynamic mapping approach for biobotic sensor networks under weak 
localization information, to obtain estimated local and global topological maps of unknown environ¬ 
ments based on only proximity information from agents who stochastically investigate an assigned 
sub-domain of the environment, with few percentage of them acting as static landmarks, while 
tracking a leading agent by maintaining their connectivity to it in order to cover, explore, and map 
an area of interest. 

The topological mapping approach used to obtain local point clouds of sub-domains can be 
viewed as reconstruction of a manifold sampled from the space of 'Di x [U-i — At, ti + At]. Hence 
the accuracy of the procedure is dependent on the number of samples taken from that space, which 
depends on the number of agents, their average velocity, and the time assigned to explore a given 
sub-domain. We will explore such dependencies to come up with bounds on the values of afore¬ 
mentioned parameters exploiting the existing literature on topological inference from point cloud 
data 43,44 . Moreover, the existence and appropriate selection of static landmarks plays an im¬ 


portant role in the construction of a distance metric under weak localization conditions. We are 
investigating the dependencies of the location and number of such landmarks on the estimation er¬ 
rors of pairwise distances, and aim to come up with criteria to choose landmarks in a more optimal 
and efficient manner. Furthermore, we are currently working on the integration of our topologi¬ 
cal mapping approach with network SLAM algorithms with range measurements to design a 
hierarchical mapping approach for disaster scenarios which require both ground and underground 
surveillance. Our future work will also include experimental results of the proposed approach with 
a swarm robotic platform, as well as the biobotic insects testbed in order to realize exploration and 
mapping with cyborg insect networks. 
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Figure 6: Simulation results for environment scenario 1 
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